Probing dipolar effects with condensate shape oscillation 
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We discuss the low energy shape oscillations of a magnetic trapped atomic condensate including 
the spin dipole interaction. When the nominal isotropic s-wave interaction strength becomes tunable 
through a Feshbach resonance (e.g. as for 85 Rb atoms), anisotropic dipolar effects are shown to be 
detectable under current experimental conditions [E. A. Donley et al, Nature 412, 295 (2001)]. 



Collective excitations play an important role in prob- 
ing microscopic interactions 0. The recently available 
gaseous atomic Bose- Einstein condensates (BEC), have 
proven to be a profitable testing ground for such studies 
0. Atomic BEC are dilute with properties completely 
determined by binary interactions. Experimentally they 
are created at very low temperatures when the short 
range atom-atom (collision) interaction can be described 
by a single parameter: the s-wave scattering length a sc . 
All higher order partial wave collisions are suppressed at 
zero collision energy in a short ranged potential. This 
enable atoms to be modelled as hard spheres of radius 
a sc , reflecting the isotropic ground state interaction. In- 
side an electric or magnetic field, however, ground state 
atoms may be polarized, e.g. in a magnetic trap the 
direction of atomic spins (of the valance electron for al- 
kali) becomes aligned. The resulting dipole interaction 
between condensed atom pairs many not simply average 
out. In this article, we investigate such spin dipole effects 
on shape oscillation frequencies of an atomic condensate. 
Although dipolar effects are typically small compared to 
the dominant s-wave contact interaction, our study in- 
dicates that these shifts become observable in the 85 Rb 
BEC H when a Feshbach resonance is used to tune 
a sc near zero ||. 

Dipolar interaction in atomic BEC leads to physics be- 
yond the usual s-wave contact term, mainly because of 
the modified low energy collision threshold behavior due 
to the anisotropic nature of the interaction S . For mag- 
netic spin dipoles, the net effect in the dilute gas sample 
is rather small, one can therefor approximate the com- 
plete two-body interaction by fg, M 



FIG. 1: Geometry for the interaction of two aligned dipoles. 



-¥ 



(a) 



(b) 



(c) 



V(R) = g S(R) + g 2 



1-3 cos 2 R 



(1) 



where R = r — r', and go = 4:nh a sc /M is the s-wave 
contact term. More generally, the dipole strength gi 
equals cv 2 (0)£ 2 (a(0) atomic polarizability) or /i 2 (// mag- 
netic dipole moment) respectively for electric or magnetic 
dipoles. Figure [l] illustrates the geometry for two aligned 
dipoles along the local magnetic field (z-axis). Several 
interesting properties of a dipolar condensate have al- 
ready been discussed ||, 0, [|, gj. In this paper, we focus 
on the shape oscillations of a trapped dipolar condensate 
assuming a tunable a sc . 

In the standard approach, low energy collective exci- 
tations are described by the Bogoliubov theory |1Q] . The 



FIG. 2: Collective excitations of a condensate with cylindrical 
symmetry. Mode a is purely radial due to the cylindrical 
symmetry. It's angular momentum projection along the z- 
axis m — 2. Mode b and c are respectively the quadrupole 
and monopole oscillations. They are customarily called the 
low and high m = modes. 



condensate wavefunction <fi(r, t) is governed by the Gross- 
Pitaevskii equation, while the non-condensed atoms are 
described by quasi-particles |ll|| , some of which have been 
studied experimentally [ p"2[ |13| . The inclusion of the non- 
local dipolar interaction makes the Bogoliubov approach 
impractical to implement numerically. We therefore will 
rely on two alternative methods to study the three char- 
acteristic shape modes |l2| as graphically illustrated for 
a cylindrical symmetric trap in Fig. ^. 

First we adopt the highly successful time-dependent 
variation approach used in Ref. jw| by assuming a Gaus- 
sian ansatz 



<t>(x,y,z,t)=A(t) J] 



= -»? 2 /2^+i») 2 /3,(t) 



(2) 



r]—x,y,z 



For an harmonic trap V t = MJ2 v =x,y,z hl ? l 'ri 2 /2 { v v = 
wAjj), the equations for variational parameters are 
equivalent to the classical motion of a particle (with co- 
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FIG. 3: Mode frequencies for N = 10 4 and A = 1. Solid 
(dashed) lines are results excluding (including) dipole inter- 
actions, ao stands for Bohr radius. 



ordinate q^) inside an effective potential 



U(q x ,q y ,q z ) 

ri=x,y,z 
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(3) 



where N is the number of condensate atoms. The equilib- 
rium location (g°) of Eq. (|J) then determines condensate 
size, while linearized shape oscillation frequencies are de- 
termined by the second order derivative U vri i(q x ,q v ,q z ) 
evaluated at q®, For a cylindrically symmetric trap (A x = 
A. y = 1), we take q x = q y = q r and A 2 = A. The last 
integral in (H) as well as U vv ' all become analytically 
computable H. The resulting U nn ' matrix is symmetric 
(Urm' — U ri ' rl '), diagonalization of which gives the three 
mode frequencies v a = ■sfU\\ — V\i and 
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± 



u 2 
^33 



) 2 + 



8U 2 3 



1/2 



(4) 



In the JILA experiments ]|, || with 85 Rb in the \F — 
2, Mp = 2) state, the valance electron spin gives rise to 
an aligned magnetic dipole moment of /1 = 2/^^/3 ([Ab 
being the Bohr magneton) . We take the radial frequency 
to be v r — 17.35 (Hz), in addition to a tunable a sc and 
trap aspect ratio A as in the experiment H. Our re- 
sults confirm that spin dipole effects become detectable 
in terms of shifts of the shape oscillation frequencies. We 
first report variational calculated results as the analytic 
formulae obtained allow for a careful analysis of the un- 
derline physics. 



Figure H shows the three mode frequencies for N = 10 4 
in a spherical trap (A = 1). The solid and dashed lines 
denote respectively mode frequencies without and with 
the dipolar interaction. The vertical lines show the criti- 
cal values of a£ when the mode character switching b <-> c 
occurs. This mode switching occurs whenever the overall 
mean-field condensate interaction changes from repulsive 
to attractive or vice versa. When squeezed along the ra- 
dial direction, mode c (b) became predominantly excited 
for an overall attractive (repulsive) condensate as atoms 
are pulled in (pushed out) along the z-direction. With- 
out the spin dipole interaction, this mode switch always 
occurs at a sc = 0. The anisotropic dipole interaction af- 
fects the overall condensate mean field and the stability 
depending on the trap aspect ratio A (|, ||. Depending 
on the configuration of the dipoles, dipole-dipole inter- 
actions can be either attractive or repulsive. For two 
dipoles, if they were placed in a plane perpendicular to 
their polarization (TT); they repel each other. On the 
other hand, they attract each other if they were placed 
along the direction of their polarization (—>—>). For a 
cloud of trapped dipoles one therefore expects that re- 
pulsive interaction increases as one increases A, which 
leads to increased condensate stability. In Fig. ||, we 
display the A dependence of the critical scattering length 
a^ c as obtained from the variational calculation. We find 
that the effective scattering length a e g (introduced ear- 
lier by us in ||) provides a reasonable approximation to 
the sign of the overall condensate mean field. In terms of 
the actual values, however, the results in the figure show 
that a e ff differs significantly from a^ c when A deviates 
from the neighborhood of unity. 

In Fig. [|, we show the A-dependence of the dipole 
induced fractional changes to the mode frequencies for a 



condensate with 10 atoms at a s . 



-l(ao). For strongly 



prolate or oblate traps, the shifts are in a few percentage 
range even with such small numbers of atom. Figure |6| 
summarizes our results for the atom number dependence 
of the dipole induced frequency shifts for the trap as- 
pect ratio A = 6.8/17.35 g. It was shown earlier at 
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FIG. 5: Fractional change (%) of mode frequencies for a sc — 
— 1 (ao) and N = 10 4 . Similar results are obtained for a BC — 
and 1 (ao). 
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FIG. 6: Atom number dependence of the dipole interaction 
induced shape oscillation frequency shifts. Mode a, b, and 
c are labelled respective with solid, dashed, and dot-dashed 
lines. 



increased values of N, a dipolar condensate always col- 
lapses irrespective of the sign of a sc . The Gaussian ansatz 
(| ) becomes questionable near collapse as reflected in the 
seemly divergent results when N is increased. 

These results (Figs. |[^) clearly show that dipolar ef- 
fects are detectable in the current 85 Rb BEC setup, given 
the extraordinary capability of 0.1% frequency measure- 
ment To confirm the validity of the variational cal- 
culations, we have invested considerable effort in an exact 
numerical method based on the time dependent Gross- 
Pitaevskii equation p7|. By applying external periodic 
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FIG. 7: Typical numerical results for condensate width and 
its corresponding Fouier transformed signal. Mode b and c 
are clearly identified. The parameters are a sc = 0, A = 1, and 
N = 40000. 



forcing terms described by the potential 



rj—x,y,z 







with A„(t) =ryo+ 8„wi(Cl„t + 0„) to the condensate 
ground state, selected shape oscillations become pre- 
dominantly excited, as first numerically implemented by 
Ruprecht et al. [ pT| . 

After a selected duration T, typically several periods 
of the trap radial oscillation, the forcing term Vf is tuned 
off. The free propagation of the time dependent Gross- 
Pitaevskii equation is continued, and the dynamic con- 
densate width yj (i] 2 (t)) is sampled. The shape oscilla- 
tion frequencies are then identified by taking the Fourier 
transformation of yj (r) 2 (t)). A typical result from this 
calculation is given in Fig. [?], which shows remarkably 
clear signal. By varying V rn 770, S v , and our results 
are self-consistently checked, i.e. to be independent of 
all parameters involved as they should be in the small 
amplitude oscillation limit and to be numerically accu- 
rate. For the experimental parameter of v r — 17.35 (Hz) 
and v z = 6.8 (Hz). We have computed the shape os- 
cillation frequencies numerically for a sc — 0,1,-2 (ao), 
and compared them with the variational results in Table 
[j]. The quoted errors for the numerical results shown in 
the tables mainly come from the limited sampling win- 
dow in time. Further improvement is hard as the cal- 
culations are time consuming and numerical accuracies 
become difficult to control at long times. We find that 
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TABLE I: Mode frequencies (in unit of v r ) with dipolar in- 
teraction 



a S c(ao) 


mode 


numerical 


variational 


-2 


b 

c 


2.0010 ± 0.0040 
0.8578 ± 0.0085 


2.0016 
0.8657 





b 

c 


1.9959 ± 0.0028 
0.7874 ± 0.0022 


1.9964 
0.7895 


+1 


b 

c 


0.7705 ± 0.0060 
1.9974 ± 0.0032 


0.7675 
1.9974 



TABLE II: Mode frequencies (in unit of v r ) without dipolar 
interaction 



a BC (a ) 


mode 


numerical 


variational 


-2 


b 

c 


2.0073 ± 0.0025 
0.8546 ± 0.0051 


2.0075 
0.8591 


+ 1 


b 

c 


0.7664 ± 0.0054 
2.0004 ± 0.0024 


0.7622 
2.0004 



the variational results are consistent with the exact nu- 
merical results. In order to affirm such mode frequency 
shifts are indeed from the dipolar interaction, rather than 
a mis-calibration of a sc , we also need to validate the vari- 
ational approach in the absence of the dipole interaction. 
This was explored earlier in Ref. where the numer- 
ical and variational results of the condensate dynamics 
were compared. Our results are presented in Table [Tl| It 
convincingly proves that dipolar efforts reported in the 
Table [l| is due to the physics of spin dipole interaction. 

In conclusion, we have studied low energy shape oscil- 
lations of a trapped dipolar condensate. Using a mag- 
netic field dependent Feshbach resonance to tune the s- 
wave scattering length to around zero, we have shown 



that the weak spin dipole interaction becomes detectable 
as shifts to shape oscillation frequencies under currently 
available experimental conditions || 0]. These shifts 
grow with the number of trapped atoms, and are around 
1% level with as few as 10 4 atoms. We have also inde- 
pendently verified the accuracy of the variational calcu- 
lation by developing a rigorous numerical approach for 
the three predominant shape oscillation modes. Near a 
Feshbach resonance, significant atom loss might occur as 
found in the Na experiment ]lq| . The subsequent damp- 
ing could lead to a broadening of the shape oscillation 
thus masking the direct observation of the proposed spin 
dipole effects. Fortunately for 85 Rb atoms, the condition 
of a sc (B) = corresponds to far off resonance on the 
high B-field side, where impressive controls have been 
demonstrated without any significant loss Over the 
last few years, mean field theory has proven to be re- 
markably successful when applied to BEC physics. The 
s-wave go contact pseudo-potential has made the concept 
of scattering length a sc widely popular. It is often said 
that scattering length is the only relevant atomic param- 
eter since the net interaction effect scales as -/Va sc /aho for 
a harmonically trapped atomic condensate, where ai 10 is 
the trap size. Dipole interactions as discussed in this 
article points to interesting physics beyond the s-wave. 
Within the current experimental regime of dilute atomic 
gases, mean field theory remains applicable and allows for 
the calculation of dipolar induced shifts of collective ex- 
citation frequencies. Successful experimental verification 
of our predictions will shed new light on atomic BEC. 
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